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Cosmic Ray ^He Measurements 
R. A. Mewaldt 

California Institute of Technology 
Pasadena, California 91 125 USA 

Cosmic ray J He/ 4 He observations, including a new measurement at —65 
MeV/nucleon from ISEE-3, are compared with interstellar propagation and solar 
modulation models in an effort to understand the origin of cosmic ray He nuclei. 

1. Introduction - The rare isotopes 2 H and 3 He in cosmic rays are believed to be of secondary 
origin produced by nuclear interactions of primary 'H and 4 He with the interstellar medium. 
There has recently been renewed interest in these isotopes as a result of indications from high- 
energy antiproton, positron, and 3 He observations that the origin of some primary H and He 
nuclei may differ from that of heavier cosmic rays. In this paper we report a new observation of 
low-energy 3 He, examine previously reported 3 He/ 4 He measurements at both low and high 
energies, and compare these with calculations of the expected 3 He/ 4 He ratio at 1 AU. We find 
no evidence for an excess of low-energy 3 He such as that reported at high energies. 
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2. Obse> vations - The n *w observation reported here was made with the Caltech Heavy Isotope 
Spectrometer Telescope tHIST) on ISEE-3 (now renamed ICE) during quiet-time periods from 
8/13/78 to 12/1/78. Figure 1 shows the He isotope distribution from the two highest energy 
intervals covered by HIST. This data i ults in a 3 He/ 4 He ratio of 0.066 ±0.016 from 58 to 77 
MeV/nucleon. Some of the 4 He in this energy interval is ’anomalous’ cosmic ray (ACR) 4 He, 
which has been corrected for using a decomposition of the ACR and galactic cosmic ray (GCR) 
fluxes (ll, and their time history. The derived correction factor of 1.1 2 ±0.06 results in an 
■observed' GCR 3 He/ 4 He ratio of 0.074 ±0.0 18. 




Figure i : The distnbution of 
quiet-tune 3 He (48 to 77 
MeV/nucleon) and 4 He (41 to 
67 MeV/nucleon) observed by 
the Caltech experiment on 
ISEE-3 during late 1978. 


Figure 2: Measured and calculated 3 He/ 4 He .-auos. Spacecraft 
observations: ■ This work, 1978; □ Goddard-UNH, 1972 
[21; O Chicago, 1973-1974. Balloon data: V Rochester, 1966 
I5fc £ UNH, 1972 (6); X UMd., 1972 [71; O UNH. 1977 
[81. Geomagnetic method: [> Ioffe, 1976 (ill; yX Chicago, 
1981 [9.101. 
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Figure 2 shows our new measurement along with selected other 3 He/ 4 He observations. The 
spacecraft observations [2,3,41 re from the 1972-1978 solar minimum period and are limited to 
the 50 to 100 MeV/nucleon interval v here contamination oy ACR 4 He is minimized. Since the 
reported observations include ACR 4 lie, each has been corrected as described above. Figure 2 
also includes observations from —100 to 300 MeV/nucleon by balloon-borne instruments 
[5,6,7,81 (here referred to as "the balloon observations"). As discussed in the Appendix, we 
believe that the observations as reported (and as plotted in Figure 2) have not adequately 
corrected for 3 He produced in the atmosphere, and a proposed correction (typically —16%) is 
therefor applied in subsequent Figures. At >300 MeV/nucleon the only observations use the 
geomagnetic method, including the recent Jordan and Meyer (J&M) measurement at —6 
GeV/nucleon [9,101, and an earlier result [1 1]. These experiments also fly on balloons, but are 
not subject to the same atmospheric corrections. 

3. Interpretation of ^He/^He Observations - To interpret tlu available 3 He/ 4 He data we use 
propagation calculations by J. P. Meyer [12], who calculated interstellar spectra for 'H. 2 H, 3 He, 
and 4 He for a variety of source spectra and mean pathlengths, using the standard Heaky-box" 
propagation model. The source spectra were of the form dJ/dT ® (T + U)~ 26 , where T is 
kinetic energy per nucleon and 0<U<938 MeV/nucleon. We ca'cuiated the effects of solar 
modulation on these spectra using the solar-minimum form of the interplanetary diffusion 
coefficient from Cummings et al. [131 and numerica solutions of the Fokxer-Planck equation 
including the effects of diffusion, convection, and adiabatic deceleration. Results of these 
calculations are shown in Figure 2 for source spectra with U-500. 

By comparison with the calculations (e.g.. Figure 2) each observation determines a "leaky- 
box" escape-length (X ), as shown for U-500 spectra in Figure 3. Note that the spacecraft and 
(corrected) balloon observations all favor X =6 to 7 g/cm , only the J&M measurement at —6 
GeV/nucleon indicates X C >10 g/cm . Table 1 summarizes the mean escape-lengths obtained. 
Note that the proposed atmospheric correction (itft Appendix) lowers the mean escape-length for 
the balloon observations by —1.6 to 1.9 g/cm" (depending on the spectrum), and generally 
improves agreement with the spacecraft observations. Table 1 indicates that softer source spectra 
(e.g., U-200) lead to a somewhat greater X e at low energies. This is both a propagation effect 
(see [121), and a result of the increased solar modulation required for soft spectra. 

4. Discussion - We wish to determine whether the 3 He observations are consistent with the 
propagation/modulation models derived for heavier nuclei, or whether there is evidence for He 
nuclei with a separate origin and/oj - history. For the J&M measurement at —6 GeV/nucleon we 
find a pathlength of — 1 5 ±6 g/cm" (in agreement with their value), independent of the assumed 
source spectrum. This value is significantly greater tii^n derived from the B/C or other Z>3 
secondary/primaiy ratios, which imply X =» 5.5 g/cm" at —6 GeV/nucleon [14], Thus, if the 
J&M measurement (and its interpretation) is correct, it does imply a high-energy 3 He excess, and 
a different ongin for at least some high-energy He nuclei 


E 

Figure 3: Plot of the escape length ^ 
deterrr’ned by the observations in S 

Figure 2. The mean and uncertainty r_ 
of the spacecraft (S), corrected (BC), g 1 

and uncorrected (BU) balloon obser- * 

vations are indicated. v 

a 

a 

X 

UJ 
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Table 1 - Mean Escape-Lengths (In g/cm s ) 


Source 

Spacecraft 

Observations 

Balloon Observations 
(~ 2d0 leV/nuc) 

Spectrum 

(~ 70 MeV/nuc) 

w/o atm. corr 

itm. corr 

(Ti-soo)" 1 -* 

6.3t0.9 

8 .O 1 O .6 

6.4il.l 

(T+200)"* '* 

8.911.0 

9.7t0.7 

7. 8 H. 3 

( 1 ) 

Uncertainty Includes systematic uncertainty In 


the magnitude of the atmospheric correction. 

2 

At low energies, the required ranges from —6 to —9 g/cm (see Table 1), with the lower 
value appropriate to U—500, a spectral form consistent with most studies of the propagation and 
solar modulation of Z^3 nuclei (e.g., [15,14]). An escape length of 6 g/cm 2 agrees well wi*h 
that derived from the B/C ratio at similar energies (see, e.g., [16,14]). We conclude that lew- 
energy observations of ^e/^He are in excellent agreement with the propagation and modulation 
parameters derived for heavier nuclei. 

The above conclusion agrees with most earlier studies of low-energy H and He isotopes that 
have included both propagation and solar modulation effects (e.g., [17,4,8]), but it is in marked 
disagreement with the recent interpretation of low-energy observations by J&M [9,101. They 
suggested that balloon observations at —100 to 300 MeV/nucleon were consistent with the 
Jt -15 g/cm“ escape-length required by their own measurement. After repeating their analysis in 
detail we conclude that J&M have significantly overestimated the pathlength required by the 
balloon da'a, as a result of a combinatiomof factors, and that self-consisteni interpretations of the 
low-energy data imply A. £ ( 4 He)< 10 g/cm 2 . This conclusion is independent of the magnitude of 
the proposed atmospheric correction for the balloon data, but it is strengthened by the apparent 
need for this correction, and also by the spacecraft observations. 

Measurements of ZS»3 nuclei imply a energy-dependent escape-length that decreases with 
energy above several GeV/nucleon. Figure 4 shows the expected 3 He/*He ratio for the energy- 
dependent escape-length (^Qp(E)) of Ormes and Protheroe [14], for two spectral forms. Both 
spectra can be seen to be consistent with the low-energy 3 He/ 4 He obsrr/ations and inconsistent 
with the J&M measurement. Although the U-500 curve falls somewhat above the data, it 
should be pointed out that Xqd(E) — 9 g/cm 2 at several hundred Me //nucleon, which is also 
greater than required by the B/C ratio. Thus the marginal agreement for U-500 is most likely 
the result of an inadequacy of the energy dependence of Xqp at low energies. 
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Figure 4: Comparison of the 
observations with various energy- 
dependent propagation models 
(see text). The box indicates the 
envelope of the corrected balloon 
observations. 
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Recent observations of an excess of antiprotons and positrons at high energies have led to 
several new cosmic ray origin and propagation models in which some nuclei have traversed a 
great deal of material. Such models also produce an excess of 2 H and 3 He. As an example, 
Figure 4 shows the predicted 3 He/*He ratio for the model of Cowsik and Gaisser [181, in whi^h a 
"degraded" component of cosmic rays orginates in thick' sources surrounded by —50 g/cnC of 
material. While this model is consistent with the J&M observation, it exceeds the observed 
3 He/ 4 He ratio at low energies. By relaxing the assumption that the "normal" and "thick" sources 
have the same energy spectra and composition, it might be possible to fit both the high energy 
data (including the ar tiproton data) and the low-energy 2 H and 3 He observations. 

Acknowledgements: I thank Dr. J. D. Spalding for performing most of the HIST 3 He data 
analysis, Drs. M. Garcia-Munoz, C. J. Waddington, W. R. Webber, and M. E. Wiedenbeck for 
discussions of 3 He measurements, and Dr. E. C. Stone for helpfui comments. This work was 
supported by NASA under grant NGR 05-002-160 and contracts NAS5-28441 and NAS5-28449. 

Appendix - Atmospheric secondaries are an important source of background for balloon-borne 
2 H and 3 He observations. One such contribution, which arises from the breakup of atmospheric 
N and O nuclei, leads to steeply falling spectra of 2 H and 3 He that are most significant below 
— 100 MeV/nucleon (see, e.g., [19,201). While this source has been taken into account in most 
previous studies, an additional source, due to the breakup of primary 4 He and heavier nuclei, has 
generally been ignored. Although we are not aware of appropriate cross section measurements 
for 4 He breaking up into 3 He in collisions with CNO, with 4 He + p cross sections [121 at — lOO 
to 300 MeV/nucleon —0.5 3 He and —0.4 2 H are produced per inelastic 4 He interaction. We 
might expect 4 He + CNO interactions to produce somewhat fewer 3 He and more 2 H than 4 He + 
p interactions [121, since CNO targets tend to fragment 4 He to a greater degree. As an estimate 
of the "fragmentation parameter" for producing 3 He from 4 He in interactions with CNO we take 
P 4 2”0-25±0.15, in which case a typical 3 He/ 4 Hc ratio of 0.1 at 0 g/cm^ will increase by —17% 
at 3 g/cm^ due to this process. Using this ecti uate (see also [2 1 j) we have re-corrected the 
balloon observations in Figure 2 to the top of the atmosphere. This presently uncertain 
correction might be measured with observations of 3 He/ 4 He vs. atmospheric depth. 
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IN1TJA1, RESULTS FROM THE CALTECH /DSRI BALLOCN- 
BORNE ISOTOPE EXPERIMENT 

S. U Schindler, A. Buffington? , E C. Christian, 

J. E. Orove, K. H. Lau Sc E. C. Slone 

California Institute of Technology 
Pasadena, California 91125 

1 L. Rasmussen Sc S. Lour sen 
D-inish Space Research Institute. Denmark 


ABSTRACT 

The Caltech/DSRI balloon-borne High Energy Isotope Spectrometer 
Telescope (HEIST) was flown successfully from Palestine, Texas on 14 
Hay. 1984. The experiment was designed to measure cosmic ray isoto- 
pic abundances from neon through iror. with incident particle energies 
from - 1.5 to 2.2 GeV/nucleon. depending on the element During - 38 
hours at float altitude. > 10 5 events were recorded with Z ^ 6 and 
incident energies >1.5 GeV/nucleon We present results from the on- 
going data analysis associated with both the pre-flight Bevalac calibra- 
tion and the flight data. 

1. Introduction. The experiment described here is a joint undertaking by Caltech and 
the Danish Space Research Institute A large-area (geometric factor ~ 0.25 m 2 sr) 
balloon-borne instrument has been developed to measure cosmic-ray isotopic 
abundances from neon through iron, with incident particle energies from ~ 1.5 to 2.2 
GeV/nucleon. depending on the element (1). The expeiiment was first flown on 14 May. 
1964. from Palestine. Texas. Prior to flight, the detector was exposed to beams of 
carbon, neon, argon, end manganese et the Berkeley Bevalac. with the latter exposure 
providing the principal calibration for the instrument. Preliminary results associated 
with the development of mapping techniques and position- determining algorithms are 
discussed, with applicction to flight data. 

2 Instrument Description The experiment employs either Cerenkov-AE-Cerenkov or 
Cerenkov-total energy techniques for isotope resolution (see reference 1). depending 
on whether the incident particle traverses the entire detector, or stops at an 
intermediate position in the instrument. Figure 1 shows a schematic diagram of the 
detector A stack of twelve Nal(Tl) disks (each nominally 2 cm thick. 52 cm diameter; 
87.2 gm/cm 2 total thickness) directly measures the energy change AE for an incident 
particle Two Cerenkov counters (Cl and C2), measure the change In Lorentz factor. 
^7 = 7. -J 9 . for the event. For a stopping particle = 1.0. Mass M is obtained from 

M = AE/A7 

Each disk making up the stack is viewed by six individually digitized 
photomultiplier tubes (PMT) This permits not only a measurement of the energy 
deposition per layer, but through intercomparison of the six PMT responses, yields the 
particle's position in ‘bat layer (2,4) Because of the large amount of material 

1 Present address Umv of Calif orma, San Diego, CASS C-011, La Jolla, CA 92093 


necessary to slow dv m the high-energ}' nuclei, typically 85% of the incident particles 
meeting the trigger requirements will undergo charge-changing fragmentation 
reactions in the detector. For the Nal stack, these events are removed through a 
comparison of the individual layer responses. Plastic scintillators (Si and SR), define 
the geometry factor, and provide additional rejection against events fragmenting 
within the Cerenkov counters 

The top Cerenkov counter consists of a mosaic of 48 aerogel radiators of 
refractive index n — 1.1, while the bottom counter employs a combination of teflon (n 
= 1.34) and Pilot 425 (n = 1.49). Results associated with the 55 Mr Bevalac calibration 
of the aerogel counter have been presented previously (3) 

In its flight configuration, the detector and associated electronics are mounted in 
an insulated pressure vi ssel. An on-board evaporative cooling system is employed to 
maintain constant detector temperature, minimizing the need for thermal corrections 
in the data. Along with housekeeping information, the 108 digitized PMT outputs are 
recorded on-board using two video recorders, and are also transmitted to the ground 
through a telemetry link for real-time analysis. 

3. Flight Details The experiment was flown from Palestine. Texas on 14 May, 1984, 
using a 17.2 x 10 e ft 3 Winzen balloon. The instrument maintained float altitude for - 38 
hours, at a typical atmospheric depth of - 5.5 gm/cm 2 During the flight, the 
geomagnetic cutoff ranged from ~ 4.5 to 5 .fi GV. The experiment was recovered with 
minimal damage to tne instrument. 

For the flight, ~ 4 25 x 10 5 events were recorded, which included > 10 J events with 
Z it 6 and kinetic energy > 1.5 GeV/nucleon. While only a fraction of the events will be 
usable for isotope analysis, the remainder are being employed for in-flight mapping of 
the detector, gam balancing, and stability checks. 

.Analysis of the flight data has shown that with minor exceptions, ail systems 
performed as designed. The evaporative cooler maintained a constant detector 
temperature of 25.0 ± 0.5“C for the majority of the flight, with a gradient across the 
Nal stack of £ 0.2°C. Examination of threshold settings and detector gains and offsets, 
have indicated an overall average stability for the flight of better than 1%. 

4. Data Analysis. One innovation in this experiment is that the Nal stack provides both 

trajectory and energy loss information for an incident particle However, this requires 
that the response of the individual disks making up the stack be accurately mapped. 
Previous analysis results have been presented based on the November. 1982 Bevalac 
calibration, which employed a beam of 55 Mn ions, at 1.75 GeV/nucleon (4). Those 
results were based on analysis of data from two central regions of the instrument. 
Our primary analysis objective has been to use the ^Mn calibration data to generate 
full-disk response maps of each PMT viewing each stack layer, and to obtain full-disk 
maps of the response ratios with various PMT combination: or position 

determination. 

The response maps for each PMT are generated using an overlay of 1 cm bins on 
each Nal disk Typically the response of a single PUTT vanes by a factor of b to 7 
across the disk The current mapping technique uses calibration data to assign an 
average response to the center of each bin For each event, a multi-point 
interpolation process is employed to obtain correction factors from the maps, with the 
final response generated from a weighted average over ail six PMT's viewing the disk. 

Figure 2 shows a histogram of the normalized energy loss (response) for layer 1. 
averaged over the entire disk. The resolution achieved of - 2.8% FWHM is only slightly 
worse than the 2.4% previously reported for central regions of the stack (4). We 
estimate that Landau fluctuations account for approximately one-half the distribution 
width Results for layers 2 and 3 exhibit similar distributions, with a FWHM of ~ 3 0% 
Additional refinements in the mapping technique are expected to improve the energy 


resolution further 

In addition to generating full-disk response maps, worn is progressing on 
developing position determining algorithms for the disks. Position resolution is 
obtained from the calibration data, by comparing the measured event position 
obtcined from Bevalac wire-chamber data, and that inferred from algorithms 
employing ratios of disk PMT responses The use of PMT response ratios effectively 
removes the energy dependence from the position determination. Currently, slx ratio 
maps are generated for each Nal disk, -sing a 1 cm spatial grid. The maps employ 
ratios of opposite sums of two and three adjacent PMT's viewing each disk. 
Intermediate positions are obtained using an mte.*pclation process. For a given event, 
the best estimate of particle position is determined . om an algorithm which 
minimizes the difference between tha mapped ratios and the measured event ratios. 

Figure 3 shows the difference distnDution for a single position coordinate (AX), as 
determined from M Mn calibration data averaged over approximately the central two- 
thirds area of layer 1 The central portion of the distribution has a FWHM of 4mm. in 
good agreement with the results previously presented, that were limited to data near 
the stack axis (4). Similar results are obtained for the orthogonal coordinate AY 
Combining both position coordinates for layer 1 produces the dilierence distribution 
shown in Figure 4. where: 

AR = (AX 2 + AY 2 )# 

From these results, the rms position resolution for layer 1 is found to be ot order 
3mm In extending the analysis to deeper stack layers, the distribution broadens to 
~ 4.4mm by layer 5 We consider these results as upper limits to the intrinsic stack 
resolution, in that coni, .buttons from wire-chamber uncertainties and multipie 
Coulomb scattering have not yet been unfolded from the data. 

5. Discussion. While we anticipate further irrprovemen in the analysis '-esuits 
described here, our results for energy resolution are already consistent with the 
design goals for ii.ot.ope resolution m the instrument. In particular, if wo extrapolate 
these results to the full 12 layers of the stack for stopping 5o Mn, the expected 
contribution to mass error from uncertainty in AE is of order 0.3%. This uncertainty 
results in an associated mass error for the instrument cf ~ 0.1T amu for Mn at the 
Bevalac energies, with a correspondingly smaller uncertainty for the lighter elements 

Prior to the detailed analysis of the ffight deta, the maps and posit'ou- 
determining algorithms developed from the ^Mn Bevalac exposure must be extended 
to the remaining stack layers. Concurrent with this effort, we are optimizing 
techniques for trajectory determination, ain* extending the position-determining 
algorithms to the outer edges of the disks 

6 Acknowledgements This work was partially supported by NASA, urder grant NCR 
05-006-160. 
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ratios, veraged over the central 
two-thirds area of the disk, for 
55 Mn ions. 
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and those obtained from layer 1 Nal 
response ratios, averaged over the 
central two-thirds area of the disk, 
for 55 Mn ions. 
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Measurements of Ye and A r F.agmrnt' \ion Cross Sections 

K. H. Lau. R. A. Mewaidt, and E. C. Stone 
California Institute of Technology 
Pasadena, California 91125 USA 

Measurements are reported of the yields of individual isotopes of Cr to 
Co (Z = 24 to 27) resulting from the fragmentation of ^Fe, and the iso- 
tope* of Mg to K (Z = 12 to 19) resulting from the fragmentation of 40 Ar. 

1. Introduction - Recent advances in the resolution and illecting ptwer of cosmic 
ray instrumentation, have led to dramatic improvemt.its in the precision of 
cosmic ray composition measurements, both elemental and isotopic. The 
interpretation of these measurements is presently limi ted by uncertainties in the 
fragmentation cross-sections needed to correct for nuclear interactions with the 
interstellar gas. Cosmic ray propagation codes now rely mainly on semi-empirical 
cross-section formulae developed by Silberberg and Tsao (SAT), which have a 
typical uncertainty of -25% [1], 

We report here relative isotope yields from the fragmentation of -380 
MeV/nucleon “Fe and -210 MeV/nucleor *°Ar in CHa targets, observed during the 
calibration of two cosmic ray spectrometers at the Lawrence Berkeley Laboratory 
Bevalac, and compere these with calculated yields based on the S&T cross-section 
formulae [lj. Preliminary results from the 40 Ar study were reported by Lau, 
Mewaidt. and Wiedenbeck (LMW) [2]. 

2. Experimental Setup - The experimental data were obtained during Bevalac runs 
in April, 1978 (“Fe) and April, 1981 (^Ar). A description of the experimental setup 
for the ^Ar calibration is given in reference [2]; with additional details in [3], The 
setup for the “Fe run was essentially identical except that the fragmentation 
products were measured during the calibration of the CaltecI: Heavy Ion 
Spectrometer Telescope (HIST) [4], launched on 1SEE-3 In the “Fe run a 587 
MeV/nucleon beam was incident on a 5.26 g/cmf thick CH 2 target, such that the 
energy of the interactions ranged from -100 to 580 MeV/nucleon, with a mean 
interaction energy of -380 MeV/nucleon. In the ^Ar runs the beam energy was 
287 MeV/nucleon, the CHg target thickness 4.1 g/ cm 2 , and the interaction energy 
ranged from -70 to 280 MeV/nucleon with a mean of -210 MeV/nucleon. We 
estimate that -70% (90%) of the analysed interactions occurred in the CHg target 
for the Fe (Ar) runs, the balance taking place in other material including A1 and 
air. 

As discussed in LMW, because the experimental setup was designed for 
calibration purposes, absolute cross-sections could not be measured. In 
particular, the detectors intercepted only those fragments emitted within -1° of 
the beam direction. On the other hand, with the excellent mass resolution 
achieved, relative fragmentation yields can be determined, which are adequate for 
many cosmic-ray studies (see, e g., [5]). 

3. Data Analysis - For the “Fe data we use the isotope identification techniques 
developed for HIST flight data (see, eg., [0]), where we have resticted the analysis 
to events stopping in the last four detectors. Figure 1 shows examples of the 
observed mass distributions, along with Gaussian fits to the data. The excellent 
mass resolution ( ff m -025 amu) permits the isotopes “Fe, “Fe, and “Fe to be 
resolved in the presence of (much more abundant) “Fe. The relative isotope 
fractions that we observe for Cr. Mn, _‘nd Fe fragments are given in Table 1. 
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Figure 1: Maas histograms of 
Cr, Mn, and Fe isotopes due 
to the breakup of 5a Fe. 
Measured at the Bevalac 
with the Caltech Heavy Iso- 
tope Spectrometer Tele- 
scope. The mass of each 
event is estimated from a 
weighted sum of two in- 
dependent measurements. 
The observed mass resolu- 
tion is ai0.25 a mu. 

m 



The analysis of the ^Ar data set has already been described [2,3]. Table 2 
summarizes the fractional yield of isotopes (within each element) for Mg to K 
fragments (12^Z^19). The 14^Zsi8 data are essentially ths samu as those in [2] 
(with minor differences due to improved statistics and refined corrections), but 
they are presented in this form for the first time The Mg and A1 data have not 
been presented before. We find excellent agreement between our isotope 


fractions and these reported by Viyogi et 
al [7], who analyzed the fragmentation 
of “‘Ar in a C target at 213 MeV/nucleon. 
This suggests that the relative isotope 
yields do not depend strongly on the 
target material. Comparing our isotope 
fractions (Table 2) with those of Viyogi et 
al. we And no evidence for a mass- 
dependent bias within an element such 
as might be introduced by our limited 
angular coverage. 

4. Comparison with Calculated 
fragmentation Yields - Ka described in 
LMW we have used a Monte Carlo 
approach developed by M E. Wiedenbeck 
to model the experimental setup and 
calculate the expected yield of each 
isotope The calculation takes into 
ac;ount both the energy and target 


Tabic 1 - “Fe FragmentaUon 

Isotope 

Fraction of Element 

Observed 

Calculated 

“Fe 

0 009: 010 

0 498: 003 

“Fe 

0 273: 009 

0.268: 002 

“Fe 

0 058: 005 

0 238: 002 

“Mn 

0 215: 007 

0.250 : 002 

“Mn 

0 280 : 008 

0.356: 003 

“Mn 

0305:003 

0.220: 002 

“Mi, 

0 148 : 008 

0 128 : 002 

®‘Mn 

0 038: 003 

0.040 : 00 1 

“Mn 

0.010:. 002 

0 008: 001 

“Cr 

0 028: 003 

0023: 001 

“Cr 

0 088 : 005 

0 129 : 002 

“Cr 

0 294 : 009 

0.375: 002 

9l Cr 

0 315: 009 

0 280: 002 

“Cr 

0 219: 008 

0 157 ±002 

“Cr 

0 058: 004 

0 036 1 00 1 

. , i 


v 






a 
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dependence of the cross- 
sections [8] for the various 
materials traversed. Tables 
1 and U include the calcu- 
lated isotope fractions for 
comparison with the obser- 
vations 

Figure 2 shows a com- 
parison of the measured and 
calculated isotope yields 
from M Fe (normalized to the 
sum of Cr and Mn). While 
there is general agreement 
on the shape of the distribu- 
tions, there are also 
significant differences. For 
example, the calculation 
predicts Mn/Cr=0.82, while 
we find Mn/Cr=1.23±0.04, a 
discrepancy also noted by 
otheis [9,10,11] who meas- 
ured at a variety of energies 
with C, CHg, and H targets. 
Of related astrophysical in- 
terest is the result that the 
radioactive isotope ®*Mn 
fraction is lower than calcu- 
lated. 
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Figure 2: A comparison of measured 
and calculated yields, normalized to 
the sum of Mn and Cr. 


Table 2 - 40 Ar Fragmentation 


Isotope 

FracUon of Element 

Obserred 

Calculated 

«°K 

0 415* 080 

0 392* 026 

*K 

0.585* 080 

0 409 * 026 

“k 

<0.03 

0 199 * 021 

»Ar 

0 509* 017 

0 477 * 008 

“Ar 

0.324*. 016 

0.357* 006 

57 Ar 

0 125* Oil 

0.143 * 004 

“A.- 

0.042* 007 

0.023*002 

“a 

0.139* 013 

0 130* 004 

“a 

0.148*. 012 

0 217* 005 

"Cl 

0.296* 015 

0 248 * 005 

“a 

0253*015 

0.273* 005 

*C1 

0.115*011 

0 112*. 004 

*c. 

0 021*005 

0 020* 002 

»s 

0009* 003 

0 003*. 000 

"s 

0 042 * 006 

0 027 * 001 

*s 

0 126*010 

0.150* 001 

■s 

0 218* 012 

0 274* 001 

**s 

0 359*014 

0 377 * 002 

*s 

0 189* 012 

0 130*. 001 

*s 

0057* 007 

0 039 * 001 


Tibia 2 (cm tin lied) 


Isotope 

Fraction of Element 

Observed 

Calculated 

sap 

0 027* 006 

0.019* 001 

S4p 

0 088 * 011 

0 088*. 001 

»p 

0.241* 017 

0 300* 002 

*P 

0.341*019 

0328*002 

Sip 

0253* 017 

0.211*002 

SOp 

0.050 * 009 

0.054* 001 

“Si 

0 052* 010 

0 088* 001 

S, st 

0.137*015 

0 203* 002 

“Si 

0 387* 022 

0 408 * 002 

*Si 

0294*020 

0 206 * 002 

“Si 

0.130* 015 

0 095 * 001 

*°Al 

0040* 013 

0 052 * 001 

»A1 

0.158* 024 

0.225*003 

“A1 

0.275* 029 

0.313*003 

27 A1 

0 410* 032 

0 303* 003 

“A1 

0.117* 021 

0.107* 002 

"M a 

0 094 * 024 

0 143* 002 

“Mg 

0 306* 038 

0 388* 003 

“Mg 

0 406 * 041 

0 285 * 003 

24 Mg 

0.194* 033 

0 182* 003 





The two largest discrepancies in either data set are at oa Fe and ae K. both of 
which are ~6 times less abundant (fraction of elementl them calculated. Since 
both nuclei have 1 less neutron than a "magic number" (n=?0-l for M K. n=28-l for 
^e), their yield may be suppressed if neutron emission is involved. The S*T 
formulas do not take nuclear shell-structure into account except for a "pairing" 
correction. Although there are other nuclei with n=19 or n=?7 that do not exhibit 
such a dramatic effect, we suggest that nuclear shell structure should be 
examined carefully in any attempts to improve cemi-empirical cross-sections 

Other significant discrepancies between the observed and calculated 
fractional yields (e.g., a 'Cl “Cl, ’•Cl, and M S) involve peripheral reactions, which 
S&T calculate with special formulae. The agreement for such reactions is better 
for “Fe, but these reactions should be less important at the higher energies 
appropriate to the “Fe data set. 

For both ^Fe and ^Ar fragments the medians of the observed mass 
distribute. is are lower than calculated (Figure 2 above and Figure 3 in [2]). There 
are also differences in the widths of the distributions that are less easily 
characterized. For both data sets, the ratio of the calculated to measured isotope 
fractions exhibit rms differences of 25%. Althc th consistent with the claimed 
accuracy of the S&T formulae, this demonstrates the need for further cross 
section measurements if the potential of cosmic ray composit:’ 'n measurements is 
to be realized. 

Acknowledgements: We are grateful to Dr M. E. Wiedenbeck, who developed the 
Uonte Carlo approach for evaluating the senh- empirical cross-section formulae, 
for helpful discussions on a number of aspects of this work. Dr. J. D. Spalding was 
responsible for developing most of the techniques for resolving Fe isotopes in HIST 
and offered advice on several occasions. We also thank Dr H C. Crawford for help 
with the Bevalac calibrations. This work was supported in part by NASA under 
grant NCR 05-002-160 and contract NAS5-28449. 
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Hexagonal Uniformly Redundant Arrays for Coded- Aperture Imaging 


M. H. Finger and T. A. Prince 
California Institute of Technology, Pasadena. CA 91125 


1. Introduction. 


Uniformly redundant arrays are used in coded-aperture imaging, a technique for forming im- 
ages without mirrors or lenses. This technique is especially itr jortant for the high energy x-ray 
and y-ray region above 20 *c V. ui this technique, a mask consisting of ooaque (closed) and tran- 
sparent (open) areas is placed between the photon sources to be imaged and a position sensitive 
detector or a detector array. Each source casts a shadow pattern of the mask or aperture onto the 
detector. This shadow pattern may be viewed as an encoded signal for that source direction. If 
each possible source code is unique, the detected composite of overlapping shadow patterns may 
be decoded to produce an image of the source distribution. 


Figure 1 shows a mask suitable for imaging 
This mask consists of an array of open (white) 
and closed (gray) cells arranged in a periodic 
pattern. The unit pattern is outlined. The mask 
in figure 1 is a uniformly redundant array 1,2 
(URA). URAs have an especially desirable pro- 
perty: the overlap between two source codes is 
independent of the source directions as long as 
the sources are sufficiently separated. Except for 
periodicity, this guarantees a unique decoding of 
the composite shadow pattern with a maximal 
immunity to statistical noise. 

To date, most work on URAs has concen- 
trated on those constructed on rectangular lat- 
hees. In this paper we focus on URAs construct- 
ed on hexagonal lathees, although many of the 
results are independent of the lathee type. 

We will present complete details for the con- 
struction of a special class of URAs, the skew- 
Hadamard UfLis, which have the following pro- 
perties: 



Figure 1. An HURA of order 79. 


1 ) They are nearly half open and half closed. 

2) They are antisymmetric (exchanging open and closed cells) upon rotation by 180* except 
for the central cell and its repetitions. 

Some of the skew- Hadamard URAs constructed on a hexagonal lattice have additional sym- 
metries. These special URAs that have a hexagonal unit pattern, and are antisymmetric upon ro- 
tation by 60*. we call hexagonal uniformly redundant arrays (HURAs). The mask in figure 1 is 
an HURA. 


HURAs are particularly suited to our application, y-ray imaging in high background situa- 
nons. In a high background situation the best sensitivity is obtained with a half open and half 
closed mask. Furthermore, systematic variations of the detector background from position *o po- 
sition can be larger than the variations in detected flux due to sources. With a skew-Hacumard 
URA a simple ro*ation turns the mask into a near anti-mask, allowing exact position-by-position 
background subtraction. Also, the hexagonal symmetry of an HURA is more appropriate for a 
round position-sensitive detector or a dose-packed array of detectors than a rectangular sym- 
metry. This is especially true for shielded detector systems where compactness is at a premium. 
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2. Mathematical Structure of URAs 

A URA is defined within a unit pattern which is repeated periodically. The number of cells 
in this unit pattern is the order v of the URA. Of these cells, k of them are closed and v-k are 
open. The uniform redundancy property of URAs involves how fr. quently a given displacement 
between closed cells occurs. We will consider a cell within a repetition of the unit pattern as 
equivalent to the corresponding cell in the unit pattern, and will therefore define the difference 
bet-tet *—o cell centers as the vector displacement between th a, folded back into the unit pat- 
tern. For a URA, all possible differences occur a uniform number k bmes among the pairs of 
dreed ~'l centers. This property guarantees the uniform overlap of Lource codes discussed in 
the introduction 1 . 

The mathematical structure of a URA is that of an Abelian group difference set 3 , which is 
specified by an Abelian (additive) group G of order v, and a set D of k elements of G with the 
property that any possible nonzero difference occurs exactly k times between elements of D. For 
a URA the group G is the lattice translations modulo the periods of the mask pattern, and the set 
l> contains those translations that take the central cell to a closed cell. The simplest examples of 
group difference sets are one-dimensional sets known as cyclic difference sets defined on the 
group of integers mod v. These play an important role in the construction of many URAs. 

URAs in the class consided in this paper, the skew-Hadamard URAs , are nearly antisym- 
metric. That is, for any nonzero element in the group G, either it or its negative but not both, are 
contained in the difference set D 3 . These skew-Hadamard URAs are a subset of the Hadamard 
URAs which are nearly half open and half closed. Hadamard URAs are characterized by the 
parameters v-4n-l, lc-2/i-l, X-n- 1 for some integer n. 

Johnsen 4 has proven two interesting facts about skew-Hadamard URAs : 

1 ) All skew-Hadamard URAs have a cyclic group G, and therefore can be constructed from 
skew-Hadamard cyclic difference sets. 

2) All skew-Hadamard cyclic difference sets are of prime order v - 3 mod 4 and can be gen- 
erated from the quadratic residues mod v . 

These facts allow us to present a construction for all antisymmetric or skew-Hadamard URAs. 

3. Construction of Skew-Hadamard URAs from Quadratic Residues 

We now pres e nt a simple procedure for gen- 
erating any skew-Hadamard URA. Ac example 
constructed on a hexagonal lattice is shown in 
figure 2. The procedure consists of the following 
steps: 

1) Choose the lattice on which the URA is 
to be constructed. The lattice is defined 
by picking two basis vectors ? 0 and . 

For our example we have chosen a hex- 
agonal lamce, which has the basis vec- 
tors separated by 60* . 

2) Choose as the order of the URA a prime 
of the form v - 4n-l . In our example 
v - 23 . 

3) Construct the order v skew-Hadamard 
cyclic difference set from the formula 

D-{l 2 .2 2 (-iy-) 2 ) mod v (1) 

4) Choose an integer r and label all the cells so that the cell centered at /?o+;7i is labeled 
with 

/ - (/ + rj) mod v 

and make all cells with labels in D closed. In our example r- 5. 



Figure 2. Construction of a skew-Hadamard 
URA of order 23. 


( 2 ) 
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The heart of this procedure is the construction of the skew-Hadamard cyclic difference set in 
step 3. For a proof that this is a difference set see Baumen 3 . Step 4 transfers the difference set 
properties onto the lattice. This is done through the labeling, which transforms addition mod v to 
vector addition on the lattice modulo the resulting lattice periods. 

The freedom available in this procedure rests in the choice of the lattice, the choice of the 
order v,and the choice of the multiplier r The lattice type will determine what symmetries can 
occur. The possible orders form a rather dense set, the hist few choices being 
v-3, 7,1 1,19,23.31,43,47,59,67.71,79, and 83. The multiplier r determines the periods of the 
URA, and hence the shape of the unit pattern. Many of the v available choices result in URAs 
that are related by the symmetries of the lattice. 


4. Hexagonal Uniformly Redundant Arrays 

Of the large number of skew-Hadamard URAs, all of which can be constructed by the pro- 
cedure in section 3, we wish to pick out those that have a hexagonal unit pattern when construct- 
ed on a hexagonal lattice. These we call hexagonal uniformly redundant arrays (HURAs) 5 . For 
an HURA each period when rotated by 60* is again a period. It can be shown from equation (2) 
that this is only possible if the multiplier r satisfies 

r 1 - r-1 mod v (3) 

This property has a simple geometric interpretation: a cell labeled / when rotated by 60* will have 
the label rl mod v. This feature, and the properties of quadratic residues modulo a prime, causes 
this restricted set of URAs to have a rotational antisymmetry upon rotation by 60* as well as 
180*. 

It can be shown that HURAs exist for order v-3 and any prune order of the form 12/1+7. If 
HURAs related by symmetry are considered equivalent, then for each of these orders there is a 
single HURA. The number of available HURAs is still large; figure 3 shows the number of 
HURAs with order v or less for v up to 10,000. In figures 4 through 8 we show examples of a 
few moderate order HURAs. 




Figure 3. The number of 
given order or less. 


HURAs with a 


Flew* 4. An HURA of order 67. 
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Figure 5. An HURA of order 139. Figure 6. An HURA of order 151. 



Figure 7. An HURA of order 331. Figure 8. An HURA of order 619. 


WhJe HURAs are mathematically interesting constructs, they also have numerous attractive 
features tor applications in astronomy instrumentation. For instance, the HURA of figure 1 is 
being implemented as a 115 kg lead coded-aperture mask on a Caltech imaging y-ray telescope 
(see OG9.2-2). 
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A Balloon-Borne Imaging Gamma-Ray Telescope 

W. E. AJthooM!, W. R. Cook. A. C. Cummings, M H. Finger, 

T. A. Prince, S. M. Schindler, C. H. Starr, and E. C. Scone 

California Institute of Technology Pasadena. CA 91125 USA 

1. Introduction. This p*per describes a balloon-borne coded -aperture 7 -rav telescope for galactic 
and extragalacuc astronomy observations. The instrument, called GRIP (Gamma Ray Imaging 
Payload), is designed for measurements in the energy range from 30 keV to 5 MeV with an angu- 
lar resolution of 0.6* over a 20* field of view. Distinguishing cnaractenstics of the telescope are a 
rotating hexagonal coded-aperture mask and a thick Nal scintillation camera. Rotating hexago- 
nal coded-apertures and the development of thick scintillation cameras are discussed in Cook et 
al. (1984 and 1985, referred to as Papers I and II respectively). 

2. Instrument Description. The basic elements of GRIP are shown in Figure 1 The telescope 
consists of a shielded detector system separated by 2.5 m from a lead coded-aperture mask. The 
pnmay detector is a position-sensitive scintillator which records the characteristic spatial pattern 
of photons cast by a y-ray source through the mask. 

The mask is trade of lead hexagons 2 cm thick and 2.5 cm across ( fLat-to-Qat). supported 
by an Al honeycomb sandwich which is transparent at y-ray energies. The aperture contains 2000 
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Table I: GRIP Balloon Telescope 


Primary 

Detector 

- - — i 

41 cm x 5 cm Nal Anger Camera 
Position Resolution < S mm rms 
(0.1-5 MeV) 

Shield 

Back Plate 5 cm Nal 
Sale 16 cm plastic xanullaior 

Mem 

Hexagonal L'RA 2000 cells (2.5 cm) 
Rotation Rate: 1 rpm 
Spacing: 2.5 m from Nal deteaor 
Size 1.2 m diameter x i cm (Pb) 

Energy Range 

C 03- 5 MeV 

1 Energy 
Reaoluuoo 

S 3 keV FWHM • 50 keV 
70 keV FWHM • 1 MeV 

Imaging 

Resoluuon 0.6* 

(1070 puds in 20* FOV) 
Angular Localization 3 arc min 
(lOo source) 


Figure I. 
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hexagonal cells of which half are open and 
half contain a lead hexagon. The cell oattern 
(see Figure 2) forms a hexagonal uniformly 
redundant array (HURA) that u opuraal for 
coded- aperture imaging. HURA’s are dis- 
cussed in more detail in Paper I and in paper 
OO 9.2-1 in these proceedings. 

Continuous mask rotation imposes an 
additional level of coding on tne y-ray signal. 

Due u> the antisymmetry of the coded- 
aperture pattern under 60 degree rotation 
(open and closed cells inter. nange for all but 
the central cell) the y-ray signal at each posi- 
tion on the detector is ume-modulaied with a 
50% duty cycle. This feature allows a com- 
plete background subtraction to be per- 
formed for each detector posiuon once every 
20 seconds assuming a I rpra rotation rate. 

In addition, the continuous rotation permits 
extension of the held of view to 20 de gr ee s , 
greatly increasing the number of pixels im- 
aged [Paper I). 

The primary detector is a Nai(Tl) 
came, a plate 41 cm in diameter and 5 cm 
thick manufactured by the Harshaw Chem.'oJ Co. The NaJ is viewed by nineteen 3 inch 
Hamamatsu R1307 photomultiplier tubes (PMT ■ ) which are individually pulse height anafzed. 
The PMT gains are calibrated continuously using pulsed LED’s for snort term relauve gain cali- 
bration and an ’ 4l Am tagged y-ray source for long term absolute gain calibration. The Am 
source is situated I m above the coded-? perrure mask and can be imaged continuously during 
flight, allowing a thorough checkout of the mask -detector imaging system. 

Background suppression is provided by an anu-coincidence shield. On the side are 12 
plastic scintillator modules which form a cylinder — 16cm thick. Each module is viev-al by a 
angle 5 inch Hamamatsu R14I6 PMT The 'ower shield section is a NaJ camera plate identical 
to the primary sensor. Further background suppression is provided by the pi unary scintillation 
camera itself. The PMT pulse heights contain informauon on the depth of the interaction in the 
crystal. Tht«j the lower half of the detector can be used as an effective "integral shield" for the 
reduction of background at low energies [Paper 111. 

The telescope is mounted on an elevation pointing platform suspended from an azimuthal 
torquing system. Azimuthal stabilization and orientation are achieved using active magnetome- 
ter feedback to the azimuthal torque motor. Elevauon orientation is under command control. 
Two Schonstedt MND-5C-25 magnetometers provide 2-axis aspect information accurate to 1.5 
arc minutes. This aspect informauon is recorded in the telemetry stream and allows ccrrecuon 
of the * •'em pouuons for pointing inaccuracies such as displacement and rotarion of the telescope 
held of vim. 

For me initial flight of GRIP, we will record all nineteen 12 bit PMT pulse heights for 
each event. Event rates of up to 5x lO J /s are possible and consequently a data recording system 
with a ! Mbit/s data rate is required. We have developed a 1.4 Mb/s recording system with z5 
Gbyte capanty using commercial VCR's and audio digitizers. This development is described in 
paper OG 9.3- 1 1 in tnese proceedings 

The chararv nr'cs of the instrument are summarized ic Table I. 

3. Instrument Performance. 

Position and Energy Resolution: The GRIP scintillauon camera has been designed to have — I cm 
rWHM or better posiuon resolution over an energy range from 100 keV to 5 MeV Figures 3a 
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and 3b show histograms of y-ray event positions computed by a maximum likelihood method for 
beams of photons of 122 keV ( 57 Co) anH 662 keV ( 137 Cs) incident on the center of the detector 
[Paper III. At the lower photon energy, the 10.5 mm FWHM of the distribution is dominated by 
photon statistics. At the higher photon energy, both Compton scattering and photon statistics 
contribute to the 7.0 mm width of the distribution. The effect of Compton scattering is most no- 
ticeable in the extended uon-Gaussian tails of the distribution. 

Although the standard deviation of the distribution of computed event positions increases 
with energy, the FWHM of the distribution continues to narrow due to an increase in the yield of 
optical photons per y-ray event (Paper III. As a consequence, the point-source angular resolution 
of the GRIP telescope improves with energy, lue primary effect of Compton scattering is a 
reduction in sensitivity due to a removal of events from *he core to the tail of the position distri- 
bution. 

The energy resolution of the GRIP scintillation camera is comparable to that of single 
PMT Nal detectors. We have measured a resolution of 7% FWHM at 662 keV. 

Imaging: Figures 4a and 4b are laboratory images of 122 keV ard 662 keV y-ray sources taken 
with the fully configured GRIP telescope. The sources were suspended 10 m from the coded- 
aperture mask and the imaging algorithms were adjusted to account for the finite distance to the 
source. The images demonstrate the ability of the telescope to locate and resolve point sources at 
y-ray energies. 




DETECTOR POSITION (cm) 


Figure 3a. Histogram of computed event po 
smons showing posiuon lesoluuon at 121 
keV. 


Figure 3b. Histogram showing position reso- 
lution at 662 keV. 


Effective Area and Sensitivity: The geometrical imaging area of the GRIP scintillation camera is 
approximately 615cm 2 . This area is determined by the maximum radius (~14 cm) fr.r which 
good position resolution can be maintained. Additional factors determining the imaging effective 
area are the point-spread position determination funcuon shown in Figure 3. the full energy 
detection efficiency, and die mask contrast [Paper 11. Figure 5 shows a plot of imaging effective 
area versus energy for the GRIP telescope. 

The instrument sensitivity depends on the observed oackground which depends on such 
factors as flight location, zenith pointing angle^ shielding, and instrument mass. We estimate our 
sensitivity to be approximately 1 x 1 0* 5 ph/cnF s kev at 100 keV and IxlO* 6 ph/cra 2 s kev at 1 
MeV for a 3o 8 hour observation from equatorial latitudes 

4. Flight Plans. The GRIP telescope is scheduled for an initial flight from Palestine. Texas in Fall 
1985. Observing targets include the Cygnus region. NGC4151. and the Crab region Future 
flights are auticipated from both the northern and southern hemispheres. 
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Balloon-Borne Video Cassette Recorders For Digital Data Storage 

W.E. Aithouse and W.R. Cook 

California Institute of Technology. Pasadena. CA. 91125. U SA. 


1. introduction 

A high-speed, high-capacity digital data storage system has been developed for a new 
balloon-borne y-ray telescope (see paper OG9.2-2). The system incorporates sophisticated, yet easy 
to use and economical consumer products: the portable video cassette recorder (VCR; <ma - itla- 
tively newer item - the 'digital audio processor". The in-Qight recording system employs eight 
VCRs and will provide a continuous data storage rate of 1.4 megabits/sec throughout a 40 hour 
balloon flight. Data storage capacity is 25 gigabytes and power consumption is only 10 watts. 

2. The Digital Audio Processor 

The key to rapid development of a VCR based digital data storage system was the availability 
of the digital audio processor, a consumer product intended to allow the use of a home VCR as a 
digital audio tape recorder. A simplified block diagram of such a processor, the Sony PCM-701ES, 
is shown in Figure l. In normal stereo recording operation the incoming left and right audio sig- 
nals are each digitized with 16 bit resolution at a rate of 44 kHz. The digital data are multiplexed 
onto a single serial line and input to the 'record data processor* at a rate of 1.4 megabits/sec. The 
processor adds error detecting and correcting codes, interleaves the data, and finally generates a 
video signal in standard NTSC television format suitable for recording on any home VCR. On 
playback the digital data are unscrambled, corrected for errors, and sent to 16 bit DACs which 
reconstruct the original audio waveforms. The powerful "cross interleaved" error correction 
scheme (1,21 virtually elirnnates errors due to dropouts on the VCR tape, as is required for noise- 
free audio reproduction. 


UEET AUDIO IN 


SIGHT AUDIO IN 


lE ft AuOtO OUT 


bight audio Out 



SONY PCV1-70IES 


Figure 1. Simplified block diagram of the Sony PCM-701ES Digital 
Audio Processor. 









ns shown in Figun 1 , the circuitry within the Sony PCM-70 1 ES is conveniently divided into 
two main printed circuit boards - an analog board and a digital board. For our purpose only the 
digital board is required. It has been incorporated in our flight recording system (Figure 2) and as 
pan of a record/playback computer peripheral (Figure 3). The digital board is very easy to use. It 
operates from a single 5 volt supply, is TTL compatible and provides standard 75 ohm video input 
and output. Power consumption is 3 watts. 

Ihe Sony PCM-70 1ES may be operated in either a 14 or 16 bit mode. In the 16 bit mode 
the full resolution of the 16 bit ADCs and D.ACs arc used, but fewer error correcting bits are 
encoded. In the 14 bit mode tape dropouts of up to 32 horizontal television scan lines can be per- 
fectly corrected, while in the 16 bit mod' error conecUs>u is guaranteed only for dropouts of up to 
16 lines 13). Nonetheless, for the convenience of a simple 16 bit word-aligned format we have 
used the 1 6 bit mode. As discused below, the error rate is quite low when high-quality tape is used. 

3. Application 

In our balloon flight application, y-ray event data are written into a 2K byte buffer memory 
while a second buffer memory is read into the PCM-70 1ES digital board in the format described 
below. Each buffer can hold 62 events, and buffers are switched synchronously every 11.7 msec, 
slaved to the PCM-70 1ES data clock. Thus, —5400 events per second can be recorded with no 
dead time, while the PCM-70 1ES is supplied with a continuous data stream of 1.4 megabits/sec. 
The NTSC video signal is fed to eight Sony SL2000 portable VCRs, which are operated sequen- 
tially wun a capacity of 5 hours recording time for each VCR. Data capacities are summarized in 
Table 1. 



H>-i vCRi 


K>H VC* 2 


t-P-l VCR 3 


vioeo out- 


rOH vcro 


f-O-j VCR 4 


H>H VCR 5 


-OH VCR 6~ 


HH VCR 7 


Figure 2. Block diagram illustrating the use of the Sony PCM-70 1 ES digi- 
tal board in a balloon-borne data recording system. 



Figure 3. Block diagram showing the Sony PCM-70 1ES digital board as 
part of a mass storage computer peripheral. 
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Table 1. Balloon-Borne y-ray Telescope 
Data Storage System 

Data Rate 

1 4x 10 6 bits/sec 
175,000 bvtes/sec 
5400 events/sec 

Recordi.ig Time 
(per tape) 
(total) 

5 hrs 

40 hrs (8 VCRs) 

Dau Capacity 
(per tape) 

2.5x 1 0 10 bits 
3.2x1 0 9 bytes 
9.7xl0 7 events 

(total) 

2x10“ bits 
2.5xl0 10 bytes 
7.8x10* events 

Power Consumption 
Digital board 
Portable VCR 
(total) 

3 watts 
7 watts 
10 watts 


We have measured the error rate using Sony L750 Ultra High Grade tape recorded at "(3III" 
speed on Sony SL2000 portable VCRs. In initial tests approximately 1.5 million 2K. byte data 
blocks were recorded on each of seven video cassettes using seven different VCRs (21 gigabytes 
total). Each data block contained a 32 bit start code, a 32 bit block count, approximately 1000 
psuedo-random 16 bit numbers, a 32 bit stop code and a 16 bit cyclic redundancy check (CRC) 
code. The tapes were played back on a single VCR and each data block was checked for an accu- 
rate start and stop code, CRC error, and sequential block count. Of the 10.5 million blocks 
checked —300 were missing (due to 'playback muting’, where the PCM-701ES simply replaced 
long stretches of uncorrectable data with a fixed data word), and 19 CRC errors were detected (due 
to 'playback interpolation' where the PCM-70 1 ES replaced isolated uncorrectable data words with 
the average of adjacent words). The fraction of blocks discarded due to detected errors was 
— 3xlO -5 , a negligible data loss rate for most data acquisition applications. A sample of 600.000 
'good' blocks were checked 100% for undetected errors and none were found. Thus the undetected 
block error rate was less than 2X10 -6 , acceptable for our present application. 

The VCR data systems discussed here were used to record and playback data for the forma- 
tion of the y-ray images presented in paper OG9.2-2. In this process 56 megabytes of data were 
recovered without a single detected error. 

This work was supported in part by NASA grant NGR 05-002-160. 
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The Energy Spectrum of Jovian Electrons in 
Interplanetary Space 

S. P. Christon, A. C. Cummings, and E. C. Stone 

California Institute of Technology, Pasadena, CA 91125 USA 

W. R. Webber 

University of Hew Hampshire, Durham. NH 03824 USA 

In this paper we report on the energy spectrum of electrons with energies —10 to —180 
MeV measured with the electron telescope on the Voyager 1 and 2 spacecraft [5/one et al., 1977] 
in interplanetary space from 1978 to 1983. The kinetic energy of <*1601005 is determined by dou- 
ble dE/dx measurements from the first two deteoors (D u Df) of i. .tack of eight solid state detec- 
tors and by the range of particle penetration into the remaining six deteoors (D 3 to D g) which 
are interleaved with tungsten absorbers. 

From 1978 to 1983 (radial range —2 to —12 AU) electrons of jovian origin were clearly 
observable for electrons stopping in £>3 (E^.4 MeV) and in D 4 (Ei> 8 MeV). For eleorons stop- 
ping in £>5 (Ei> 12 MeV), the jovian flux dominated the galactic electron flux for a period of ap- 
proximately one year near the encounter with Jupiter [Christon et al., 19851. Jovian electrons 
were also observed in £> 6 (E£.21 MeV) but not in £>7 (E£28 MeV). A detailed interpretation of 
the electron variations in all energy channels depends on an accurate subtraction of background 
induced by energetic protons of a few 100 MeV. This subtraction is facilitated by laboratory cali- 
bration results at several energies. Further results on the differential energy spectrum of jovian 
electrons and limits on the maximum detected energies will be reported. 

This work was supported in part by NASA under contract NAS 7-918 and grant NGR 05-002- 
160. 
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Precision Measurements of Solar Energetic Particle 
Elemental Composition 

H Branaman and E. C Stone 

California Institute of Technology 
Pasadena, California 91125 USA 


1. Introduction 

Using data from the Cosmic Ray Subsystem (CRS) aboard the Voyager 1 and 2 
spacecraft (Stone et al. 1977), we have determined solar energetic particle abun- 
dances or upper limits for all elements with 3 ■£• Z £ 30 from a combined set of 10 
solar flares during the 1977-1982 time period Statistically meaningful abundances 
have been determined for the first time for several rare elements including P, Cl, K. 
Ti and Mn, while the precision of the mean abundances for the more abundant ele- 
ments has been improved by typically a factor of ~3 over previously reported 
values. When compared to solar photosphenc spectroscopic abundances, these 
new SEP abundances more clearly exhibit the step-function dependence on first 
ionization potential previously reported by Cook et al (1979, 1984) and Meyer 
(1961, 1P85) 

2. JbaCiWationa 

Each CRS includes four Low Energy Telescopes (LET) and two High Energy 
Telescopes (HET) employing silicon solid-state detectors and covering a combined 
incident energy range for oxygen of 3.5 - 50 MeV/nucleon A scatterplot of LET 
data from the Z = 14 - 20 charge range is shown in Fig 1 Even relatively rare, ele- 
ments such as At and Ca are clearly resolved For the rarer elements (e.g., P, Cl. 
K). the abundances were determined by performing maximum-likelihood fits to the 
rare element peak and its two usually more abundant neighbors For the most 
abundant elements (C. N, 0, Ne, Na, Mg, Al, Si. S, Ar, Ca, Cr. Fe), the abundances 
could be determined in each of the flares, these values were averaged to obtain 
mean abundances using a weighting technique that takes into account both sta- 
tistical variations and real abundance variations from flare to flare 

3. Results 

The average SEP abundances relative to silicon for elements with 3 « Z c 30 
are listed m Table 1 In Fig 2, these results for the more abundant elements are 
compared to those obtained by other recent investigations It can be seen that the 
new abundances agree well with the previously determined values but have about a 
factor of three higher precision 

4. Pise union 

An ordering by first ionization potential (FTP) of SEP composition relative to 
"solar system" or photospheric composition has been noted in the past (Hovestadt 
1974, Webber 1975, McGuire et al 1979) In particular, Cook et al (1979, 1984) 
found a step-function dependence on FTP, with elements with FTP > 10 eV depleted 
by a factor of ~5 in SEPs and elements with FTP < 8 eV approximately equally abun- 
dant in SEPs and the p biosphere This behavior may be expected on the basis of 
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Ftg 1. Scatterplot of SEP data from LET in the Z = 14-20 charge range <Z> 
is the average, and AZ the difference, of two essentially independent charge deter- 
minations obtained from the LET data for each analyzed particle 


dynamical ionization models such .is that of Geiss and Bochsler (1984) for the for- 
mation of the corona Pig 3 shows that the step-function ordering by PIP is more 
clearly defined in the new data, including the elements P, Cl, K. Ti, and Mn. How- 
ever. Pig 3 also discloses statistically significant e nha ncements and depletions of 
several elements such as Na, Pe and Ni. which are discussed by Breneman and 
Stone (1985) in the context of an acceleration/propagation fractionation of SEPs 
based on the ionic charge-to-mass ratio of the species making up the SEPs 

In Pig 4. the SEP abundances a, ompared to galactic cosmic ray source 
(GCRS) abundances (Lund 19B4). The absence of a FlP-dependence in the 
SEP/GCRS ratio indicates that the GCRS may have a PIP-dependent step-function 
fractionation like that of the SEPs (Cook et al. 1979, Meyer 1981, Lund 1984) Both 
may exhibit a depletion of elements with PIP greater than ~10 eV because of the 
smaller flux of ionizing photons with energies greater than that of Ly-a (10 2 eV) 
Both C and N appear to be distinctly different m SEPs and the GCRS 


Table 1 SEP average elemental abundances (Si = 1000) 
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• Abundances for these elements are based on fewer than 5 particles and are 
highly uncertain 
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Fig 4 SEP abundances (Si = 1) relative to galactic cosmic ray source abun- 
dances (Lund 1964). plotted vs F1P 
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Solar Coronal and Photospheric Atr indances f om 
Solar Energetic Particle Measurement'' 

H Breneman and E C Stone 

California lnititute of Technology 
Pasadena. California 61125 USA 

1. Introduction 

Solar energetic particle (SEP) elemental abc idance data from the Cosmic Ray 
Subsystem (CRS) aboard the Voyager 1 and 2 spacecraft (Breneman and Stone 
1985) are used to derive unfractionated coronal and photospheric abundances for 
elements with 3 sZs 30 We find that the ionic charge-to-mass ratio (Q/M) is the 
principal organizing parameter for the fractionation of SEPs by acceleration and 
propagation processes and for flare- to- Care variability, making possible a single- 
parameter Q/M-dependent correction to the average SEP abundances to obtain 
unfractionated coronal abundances A further correction based on first ionization 
potential allows the determination of unfractionated photospheric abundances 

2. Results 

The composition of individual flares relative to the average SEP composition 
can be described by a function that is roughly monotonic in Z but variable in mag- 
nitude from flare to flare (Cook et al 1079, 1984. McGuire et al. 1979, Meyer 1901, 



Q/M Q/M 

Fig 1. Abundances relative to the mean SEP abundance for two typical flares, 
plotted vs Q/M 
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Fig 3 Mean SEP abundance relative to the photosphere (Grevesse 1964) for the 
low-FIP elements, plotted vs Q/M The best-flt power-law has a reduced x 2 of 1.0 

1985) Using the recently reported SEP ionic charge states (Luhn et al 1984), we 
find that this compositional variability exhibits a monotonic dependence on the 
Ionic cbarge-to-m&ss ratio Q/M of the particles, as illustrated m Fig 1. Such a 
dependence is not unexpected for acceleration and propagation effects The spec- 
tral indices of the elemental differential energy spectra for a given flare also tend 
to be ordered by Q/M (Fig 2) Thus we Qerive an unfractionated coronal composi- 
tion by applying a Q/M-dependent correction to the average SEP abundances. 

3l Diacuwon 

The correction factor has been determined by comparing the SEP abundances 
of elements with low Atsl ionization potential (FIP), which display no FTP- 
dependent fractionation (Breneman and Stone 1905), to the photosphenc spec- 
troscopic abundances (Grevesse 1984) (see Fig 3). The derived correction func- 
tion. which is a power-law in Q/M with a slope of 0.66 ± 0.17, may be applied to the 
SEP abundances for all elements, resulting in SEP-denved coronal abundances 
(see Table I) As Fig 4 shows, they agree well with coronal abundances obtained by 
XUV/X-ray spectroscopy, but have much higher precision and are available for 
many more elements 

The derived coronal abundances can also be corrected for the FlP-dependent 
fractionation suggested by the dynamical ionization model of Geiss end Bochsler 
(1984). In this model, the high-FIP elements such as N, 0, F, Ne. Cl and Ar are 
depleted because their ionization times are longer than the time individual atoms 
spend in a rising spicule Thus, their abundances are corrected by the depletion 
factor of oxygen in SEPs relative to the photosphere (4 03 ± 0 26); P and S are 
corrected by the depletion factor of sulfur (1.89 t 0.17), C is corrected by the 
mean of the oxygen and sulfur depletion factors, since its proper depletion factor 
is uncertain, and the low-FlP elements, which are quickly ionized, are left 
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Table 1 SEP-derived coronal and photosphenc abundances relative to silicon 
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• Based on fewer than 5 particles and highly uncertain 


unchanged The resulting SEP-derived photospheric abundances (Table I) involve 
fewer modeling parameters than spectroscopic determinations and are available 
for some elements (e g.. C. N, Ne. Ar) that cannot be observed spectroscopically 
The main differences (Fig 5) are a significantly higher abundance of Cr (and possi- 



first ionization potential (eV) 

Fig 4 SEP-der /ed coronal abundances relative to spectroscopically- derived 
coronal abundance (Veck and Parkinson 1981), plotted vs FIP 




Fig 5 SEP-der.ved photosphenc abundances relative to spectroscopically- 
derived photosphenc abundances (Grevesse 1984) The estimated uncertainty in 
the FIP fractionation correction is indicated by the vertical brackets 


bly Ca and Ti) and a C abundance that is about half of the commonly assumed solar 
abundance The SEP-denved Fe/Si ratio is in agreement with the recent photos- 
phenc value (Grevesse 1984), which is 50% larger than the meteontic value 
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LATITUDE VARIATION OF RECURRENT FLUXES 
IN THE OUTER SOLAR SYSTEM 

S. P Qiriston and E. C. Stone 
California Institute of Technology, USA 

(Confirming Abstract) 

Recurrent low-energy (£ 0.5 MeV) proton flux enhancements, reliable 
indicators of corotating plasma interaction regions, were observed on the Voyager 
1 and 2 and Pioneer 11 spacecraft in the heliographic latitude range 2 # S to 23°N 
and the heliocentric radial range 11 to 20 AU [Christon and Stone. 1965]. After a 
penod of rather high correlation between fluxes at different latitudes in early 
1963, distinct differences developed in the fluxes during an overall flux decrease 
"Tie flux intensities returned to higher levels in early 1984 and differences in both 
thi* recurrence frequency and flux intensity persisted into 1985, as Voyager 1 
traveled to 23 AU and 25°N latitude. Intercompanson of data from the three 
spacecraft indicates that the flux differences are most likely due to latitudinal, 
rather than radial or temporal, variations. 
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DIFFERENTIAL MEASUREMENT OF COSMIC- RAY GRADIENT 
WITH RESPECT TO INTERPLANETARY CURRENT SHEET 

S. P. Christ on, A. C. Cummings, and E. C Slone 
California Institute of Technology 

K W. Behannon, and L F Burlaga 
Goddard ££>ac8 flight Center 


ABSTRACT 

Simultaneous magnetic field and cherged particle measurements from the Voy- 
ager spacecraft at heliographic latitude separations from 10° to 2l e are used to 
determine the latitude gradient of the galactic cosmic ray flux with respect to the 
Interplanetary current sheet. By comparing the ratio of cosmic ray flux at Voy- 
ager 1 to that at Voyager 2 during periods when both spacecraft are first north 
and then south of the interplanetary current sheet, we find an estimate of the 
latitudinal gradient with respect to the current sheet of approximately -0 15 ± 
0.05 % 'deg under restricted interplanetary conditions 

1. Introduction. Organization of cosmic ray flux about the heliographic equator 
should result from symmetries in the interplanetary magnetic field produced by the 
sun's rotation However, if particle drifts are imf ^rtant in cosmic ray transport, the 
cosmic ray flux should be organized with respect to the interplanetary current 
sheet. This symmetrizing property of the current sheet has been demonstrated in a 
three-dimensional model of cosmic ray transport including particle drifts [Kota and 
Jokipii, 1983]. At 1 AU, gradients with respect to the current sheet G « of up to ~ 
-0 33 X/deg have been identified for cosmic ray protons with energies >100 MeV in 
various integral energy ranges [Newkirk et al., 1984, Newkirk and Fisk. 1985] 

The interplanetary current sheet is a surface 
separating regions of magnetic field pointing generally 
toward or away from the sun In its least complicated 
configuration, the near-sun current sheet is essen- 
tially a circle tilted with respect to the heliographic 
equator As the rotating sun s magnetic field is drawn 
out by the solar wind, warping of the current, sheet 
ensues, thus allowing two spacecraft separated in lati- 
tude to sample fluxes on either side of it (Figure 1) 

The near-sun current sheet extends up to ~60° north 
N and south S of the equator and produces an inter- 
planetary sector pattern that vanes between two and 
four sectors dunng the penod of this study 
[Hoeksema, 1984] 

We determine G « using spacecraft trajectory information, the ratio of cosmic 
ray flux, and observations of the current sheet (magnetic field sector boundaries) at 
Voyager 1 (Vl) and Voyager 2 (V2), separated in latitude by 10° to 21° Such a 
differential measurement has significantly smaller non- latitudinal intensity vana- 
tions than a single-point measurement such as the type employed by Newkirk et al 
and Newkirk and Fisk The cosmic ray flux at an observation point in the helio- 
sphere is j = ; b exp(G^(r-/^))exp(Ce(9-0e,)), where Gr is the radial gradient. G 9 is 
the colatitude gradient with respect to the current sheet, is the radius of the 
boundary of the heliosphere, j* is the cosmic ray flux at assumed constant, 
is the colatitude of the current sheet, and r and 0 are the radius and colatitude of 
the observation point. Heliographic colatitudes 0 from 0* to 180° correspond to 
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Fig 1 Schematic rep- 
resentation of NN and SS 
configuration 


\ 

> 



heliographic Latitudes from 90°N to 90°S. Time senes of simultaneous observations 
at two points in interplanetary space well separated In latitude allow estimates of G 9 
to be obtained for each interval in the series These individual estimates of G 9 - 
(lnT -Gy. (r,-r 2 ))/(9 1 -9 g ), where the subscripts 1 and 2 refer to Vl and V2, respec- 
tively, and r is the ratio of the flux at Vl to that at V2. are independent of vanations 
in the separation of the observation points. Statistically weighted averages of 

= E G 9i y (E 14 ** + 

where i(J) runs from 1 to n^n^), the number of observations when both 

spacecraft are north (NN) or both south (SS) of the current shee' The first charac- 
ter of the two character identifier represents the magnetic field region for Vl, the 
second for V2 Data for NS and SN configurations are also collected and are used 
independently, and in combination with the NN and SS data, to estimate Gy. 

2. Data Selection. 

This study covers the period from day 240. 1981 until day 190, 1983 while the V1-V2 
separation increased from 1.9 to 4.7 All in radial distance from the sun and from 10‘ 
to 21° In latitude Vl traveled from 11 to 17 AU and from 8°N to 20°N V1-V2 longi- 
tude separation was «16° Daily samples of counting rates of protons with energies 
>75 MeV and a median energy -1.1 GeV measured by the nearly identical High 
Energy Telescopes of the Voyager Cosmic Ray Subsystem [Stone et al , 1977], and V 
the ratio of Vl to V2 counting rate samples are plotted in Figure 2. Calendar days 
for the NN and SS data sets are highlighted in the upper panel by underscoring We 
have compared measurements taken at both spacecraft on the same day under the 
assumptions that convective effects are unimportant on the average, consistent with 
long- interval average observations [Venkatesan et al, 1983], and that local varia- 
tions in the cosmic ray intensity, after eliminating solar flare transient effects, are 
due to variations in the distance to the current sheet [Newkirk and Fisk, 1985] 


Fig 2 Daily sam- 
ples of the 
counting rates 
of >75 MeV pro- 
tons from Vl 
and V2 (Vl is 
offset for clar- 
ity) and T the 
ratio of Vl to V2 
counting rate 
samples. Uncer- 
tainties are due 
to counting 
statistics 



The counting rates of >75 MeV and -1 MeV protons were used to eliminate 
periods possibly contaminated oy solar flare particles and Forbush decreases Pro- 
ton energy spectra and elemental abundance measurements at lower energies were 
also used for this purpose Hourly averages of the interplanetary magnetic field 
azimuth and elevation angles (Figure 3) measured by the Voyager Magnetic Field 
Experiment were used to determine the less active magnetic field periods, sector 
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polarities and transition regions between sectors Time blocks of data were 
identified when VI and V2 were in regions of magnetic field pointing generally 
toward or away from the sun along the expected field direction, which is generally 
perpendicular to the radial direction in the outer solar system and parallel to the 
equatorial plane near the equator The minimum sector length used herein is 4.5 
day6, so that VI and V2 were most likely in field regions associated with plasma flows 
from the sun's polar coronal holes This produced Data Set (DS) 1. Then days with a 
poorly A .'fined transition region at either end of the magnetic sector and days with 
azimu.ual turbulence were eliminated in order to produce DS 11 Finally, days with 
elevation angles pointing away from the equatorial region and days with elevation 
angle turbulence were collected into subset DS lib When the elevation angle 
changes by a large amount in less than a day or two as in DS lib, it is likely to 
represent a dynamical effect such as the passage of transient material or perhaps 
reconnection at a sector boundary, either of which could alter the cosmic ray gra- 
dient Those quietest days remaining are subset DS Ila 


Fig 3. Magnetic 
field azimuth and 
elevation angles 
from VI and V2 
The azimuth, X. is 
0° when the field is 
directed away from 
the sun and the 
elevation, 6. is 0 C 
when the field is 
along the space- 
craft-sun line Two 
cycles of X are 
plotted. Arrows 
highlight transition 
regions inter- 

preted as current 
sheet crossings 
Those labeled "SB" 
are s imil ar to sec- 
tor boundaries at 1 
AU. 



3. Observations and Discussion. Determination of <G*> requires knowledge of the 
value of G^.. which can be estimated from NS and SN data sets, since the latitude 
gradients. which nave opposite signs in the two hemispheres should average to zero 
Averaged ever NS and SN. <G r > MSSN = 2.21 ±0.11 2/AU, consistent with radial gra 
clients reported by [Ve.Jcalesan et al . 1963] G,. can also be determined by averag- 

ing the NN and SS data sets, assuming that G e is negligible This results in <G r > NNSS 
= 2.29 ± 0.13 2/AU, consistent with < > acs • and suggesting that the best estimate 
is the average over all 152 days of data which yields <G r > - 2.25 ± 0 08 %/AU Note 
that an uncertainty of 0.1 */AU in the assumed G,. introduces an error of only 0.01 
X/deg in the derived <G„> 

<G«> and <G«> have been calculated both without restrictions on spacecraft 
posiuon (unrestricted averagi), and with restrictions on spacecraft position with 
respect to the current sheet based upon time to or from the closest sector boun- 
dary crossing (restricted average). Requiring V2 to be closer to the current sheet 
than VI for NN data and VI closer than V2 for SS data (see e g.. Figure l). should 
enhance the measured <G«> 



SH 4.2-9 



ORIGINAL P 
OF POOR QUAL I 


Table ] <G»> and <G*> (# days of observation) 


unrt 

DS 

tirirtrd overage 
("/deg) 

<G#>ss 

(7/deg) 

<G«> 

(7/deg) 

1 

0.02 ± 0 03 (52)' 

•0 02 ± 0.04 ( 22) 

0 02 ± 0 03 (74) 

n 

0 03 ± 0 04 (29) 

-0.01 ± 0 04 (19) 

0 02 * 0 03 (48) 

Ila 

-0 02 ± 0.05 (21) 

-0 06 ± 0 03 (6) 

-0 01 t 0 04 (27) 

Db 

0.14 ± 0 09 (0) 

O.Oi ± 0 05 (13) 

0 06 t 0 05 (21) 

restricted average 



DS 

<C*>Ah 

<G*>«s 

<G«> 


(r./deg) 

(7. /deg) 

(7. /deg) 

I 

-0 12 ± 0 04 (21) 

0 09 ± 0 07 (7) 

-0.11 ± 0 03 (28) 

11 

-0 18 ± 0 05 (10) 

0.12 ± 0.09 (5) 

-0.16 ± 0 04 (15) 

Ila 

-0 18 ± 0.05 (10) 

-0 06 ± 0 08 (1) 

-0.15 ± 005 yll) 

lib 

(0) 

0 20 ± 0 09 (4) 

-0.20 ± 0 09 (4) 


Table 1 lists values of <G*> and <G*> for both posi- 
tional entena (the lower case subscript is for the space- 
craft closer to the current sheet) using G,. = 2.25 %/AU 
for the radial gradient correction. Note that the unres- 
tricted <G«> display? no observable gradient, while the 
restricted <G«> is statistically significant and is approxi- 
mately -0 15 ± 0 05 7/deg Since differences between DS 
Ila and DS lib are statistically insignificant, we have no 
evidence that field elevation turbulence affected our 
measurements Figure 4 shows the separation of the Nn 
and sS subsets for DS II with the restricted average Note 
that, even with such a small data set, a separation of 
means and modes of the Nn and sS distributions, attribut- 
able to a small, negative latitude gradient, is evident. 
Further studies are required to determine whether such 
gradien.s are typical jf other time periods 



Fig 4 <G«> for DS 
D. restricted aver- 
age 
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Changes in the Energy Spectrum of Anomalous Oxygen and Helium 

During 1977-1985 

A. C. Cummings and E. C. Stone 

California Institute of Technology Pasadena CA 91125 USA 

*A. R. Webber 

( niversity of Sr* Hampshire. Durham. SH 03824 USA 

We have used data from the cosmic-ray experiment on the Voyagei spacecraft to 
measure the energy spectrum of anomalous O and He during the period 1977 to 
1985. We find that these spectra change dramatically after the middle of 1980, 
with the peak or plateau region of the differential spectrum shifting to a higher 
energy . This change appears to be related to the reversal of the solar magnetic 
fie'd and could irayly that particle dnfts arc important to the acceleration or pro- 
pagation of these particles. 

1. Introduction. The study of the anomalous component of cosmic rays over the solar cycle may 
prove to be key to the understanding of the role of dnfts in cosmic-ray modulation. Pesses et al. 

(1981 ) have suggested that the acceleration site of the anomalous component is most likely in the 
polar regions of the solar wind termination shock. In their model the particles drift lab* ally 
toward the neutral sheet from the polar regions dunng the last solar cycle, but dr radii' • in- 
ward along the neutral sheet during the current solar cycle. Consequently, they p t ong 
dependence of the intensity of the anomalous component near the sohr equator on solar magnet- 
ic field polarity. In this analysis we make use of measurements of the spectra of C and He from 
the Cosmic Ray Subsystem (CRS) LS/one et al.. 1977) on the Voyager 1 (VI) and 2 (V2) space- 
craft to address whether or not these spectra exhibit changes associated with the polarity of the 
solar magnetic field. 

2. Observations. The general features of the 1 1-vear cosmic-ray modulation cycle are evident in 
the counting rate of the particles which penetiatr the high energy telescope on VI, shown in Fig. 

1 for 1977-1985. In order to minimize the contamination by solar and interplanetary energetic 
particles, quiet times were selected by setting limits to the maximum low-energy He flux, in a 
manner similar to that described in Cummings et al [ i 984). The six quiet-time periods that 
were selected for analysis for V 1 are shown as horizontal bars in Fig. 1 and are labeled A through 
F. Also shown are the approximate times of the reversal of the solar polar magnetic nclds I Webb 
et al.. 1984). Note that the field reversal occurs approximately between periods C and D. 

The anomalous O spectrum for the Voyager measurements is derived from the observed 
total O spectrum for a particular time interval by subtracting a low-energy solar or interplanetary 
component and the high-energy galactic cosmic-ray component. This subtracuon procedure is il- 
lustrated in Fig. 2 for two extreme examples. The low-energ component is scaled and extrapo- 
lated from a power-law fit to the observed He spectrum in the energy range 3-6.1 MeV/nuc using 
a rauo He/O- 100 ± 50 (see Gloeckler et al (1979)). The galactic cosmic-ray O component is es- 
timated by normalizing the observed carbon spectrum to the O .ntensity in the 66- 125 MeV/nuc 
energy interval. 

FiguT 3 shows the spectra of anomalous O for the six quiet-time intervals A through F 
derived from the V2 observations in the manne- just described. The a and b panels show the 
spectra for the three intervals before and after the solar magnetic field reversal in 1980. respec- 
tively. The spectra for intervals A. B. and C are essentially monotonic with a tendency to flatten 
into a peak or plateau at less than —6 MeV/nuc. It is clear from Fig. 3b that the energy depen- 
dence of the spectra for the intervals after the held reversal (D, E. and F) is significantly different, 
having a peak or plateau at — 10-20 MeV/nuc. 

A similar change in the energy spectrum is also apparent in the VI data. In Fig. 4a we 
show the V 1 observed O spectra for intervals B and F. The galactic cosmic-ray intensity is dom- 
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F'g. 1. Three-day average counting rate of penetrating particles in the PET 1 (High-Energy Tele- 
scope number 1) of the VI CRS instrument from 1977 to 1986. The ratr is dominated by protons 
with energy > 75 MeV The bars and arrows are described in the text 



Fig 2. Quict-ume energy *pectra of O from Voyager data for two sample ume intervals The ob- 
served spectra aie shown as the solid circles. Estimated spectra of interplanetary and galactic com- 
ponents are shown as solid and dashed lines, respectively. The spectra of anomalous O are indicated 
by the open squares which arc joined by the dotted lines 
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Fig. 3. Spectra of anomalous 0 from V2 for the six quiet-time intervals. Panels a and b show the spectra 
for periods before and after the solar magnetic field reversal, respectively. The smooth curves are drawn to 
aide the eye and are not functional fits to the data. The dashed fine labeled ’A’ iu panel b is copied from 
panel a to facilitate comparison. Note the higher energy of the peak (or plateau) intensity of the spectra in 
panel b when compared to those of panel a. 


man? in the 50-125MeV/nuc energy range for both time intervals. The intensity in this energy 
range is only — 20% higher for interval 3 than for interval F, indicating that a similar level of 
modulation has been reached. Similarly, the anomalous C intensity in the 10-30MeV/nuc ener- 
gy range ; s approximately tne same for the two time intervals. However, below lOMeV/nuc 
there is a factor of 10 difference in intensity producing a striking difference in the spectrum. 

This change in the energy speefum is also evident in the He data. We show in Fig. 4b the 
observed He spectrum for the same two intervals as for Fig. 4a. Below ~8MeV/nuc a low- 
energy solar or interplanetary component is dominant and causes a sharp upturn in the spec- 
trum. For interval B the flat region of the spectrum from — 8to50MeV/nuc indicates that 
anomalous He is the primary component, dominaung the flux of solar and galactic cosmic rays. 
The period F spectrum has a much lower intensity than the period B spectrum beiow 
— 30MeV/nuc suggesting that the anomalous He spectrum has shifted to higher energies. If the 
particles are singly ionized, the energy dependence of ♦he He spectrum is expected to be obtained 
by shifting the anomalous O spectrum by a factor of —4 in energy [Cummings el al., 1984], so 
that the absence of anomalous He below 30MeV/nuc is corsistent with the reduced intensity of 
anomalous O below 8MeV/nuc. 

3 Discussion. If the changes in the energy spectra of the anomalous component are to be ex- 
plained by "conventional" modulation theory, in which the effects of diffusion, convection, and 
adiabatic deceleration are considered, and the effects of drifts ignored (see Fisk [1980], for a re- 
view of solar modulation theory), then it would require either a very prolonged "hysteresis" effect 
or a significant change in the rigidity dependence of the diffusion coefficient between halves of 
the solar cycle. A hysteresis effect (phase-lag between intensity variations of low-rigidity and 
high-rigidity particles) has been reported for the anomalous 0 component by Klecker ei al 
[1980] for the 1974-1975 period. They found a pnase lag of 72 days between the intensity of 
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Fig. 4. Spectra o f observed O (a) and He (b) for time intervals B and F from VI. 


anotralous O with 7.6-24 MeV/nuc and the intensity of galactic cosmic rays with energy 
> iOGeV. If the energy spectrum changes in our study are to be ascribed to hysteresis, then a 
phase lag of — 4 years would be required between particles with only mtdestly different rigidities 
and velocities. If changes in the diffusion coefficient between halves of the solar cycle are to ac- 
count for our observations then the spectra of other cosmic-ray species would show pronounced 
effects at the corresponding rigidities. 

It appears more likely, therefore, that the observed time variation of the energy spectra of 
anomalous O and He may be related to changes m the acceleration of panicles a: tne termination 
shock due to the change in the polarity of the solar magnetic field. Recently, Jokipii [1985] has 
used a model of the anomalous component which includes acceleration at the termination shock 
and particle anfts and has calculated spectra of anomalous O and He in the two halves of the 
solar cycle with different field polarity that resemble the observations reported here. 

4. Acknowledgments. We are grateful to R. E. Vogt and other CFS team members at the Cali- 
fornia Institute of Technology, the Goddard Space Flight Center, the University of Arizona, and 
the University of New Hampshire. This work was supported in pan bv NASA under contract 
NAS 7-918 and grant NGR 05-002-160. 
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SOLAR CYCLE VARIATIONS OF THE ANOMALOUS COSMIC 
RAY COMPONENT 


R. A. Mewaldt and E. C. Stone 

California Institute of Technology 
Pasadena. California 91125 USA 


The intensity of the "anomalous” cosmic ray component, consisting 
of He. N. 0, and Ne. has long been known to be especially sensitive to the 
effects of solar modulation. Following its discovery in 1972. this 
component dominated the quiet-time flux of cosmic ray nuclei below -30 
MeV/nucleon during the 1972-1978 solar minimum, but then became 
essentially unobservable at 1 AU with the approach of solar maximum 
conditions One recent theoretical model predicts substantial 
differences in the intensity of the ano mal ous fluxes from one solar 
minimum period to the next because of the reversal of the solar 
magnetic field. Using data from the Caltech experiments oa IMP-8 and 
ICE (1SEE-3), we report on the intensity of anomalous 0 and He at 1 AU 
during the years 1972 to 1985 !n particular, we hope to determine 
whether the anomalous fluxes will return to their 1972-1978 levels, as 
predicted by spherically symmetric modulation models, or whether they 
will fail to return to 1 AU. as suggested by modulation models in which 
gradient and curvature drifts do min ate 

Our preliminary analysis of data from 1984 shows that the intensity 
of 8 to 27 MeV/nucleon 0 is still more than on order of magnitude below 
its 1972 to 197B levels, while the intensity of 25 to 43 MeV/nucleon He is 
a factor of ~8 below its maximum level in 1977 Data from 1985 is now 
being analyzed. 

This work was supported in part by NASA under grant NGR 05-002- 
160 and contracts NAS5-28441 and NAS5-2B449 
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We find that the radial gradient of anomalous O remains constant during 1977-85 
at — 10-15%/AU although the intensity changes by more than a factor of 100. 

These results can be used to deduce that most of the modulation of the intensities 
of these particles is occurring beyond 27 AU. We also find evidence for a latitudi- 
nal gradient of — +3%/degr?e at low energies (7.1- 10.6 MeV/nuc). 

1. Introduction. In this analysis we make use of measurements of the spectrum of O nuclei from 
the Cosmic Ray Subsystem (CRS) on the Voyager 1 (VI) and 2 (V2) spacecraft [Stone et al.. 
1977] and from the Goddard-University of New Hampshire experiment on Pioneer 10 
[ McDonald et al.. 1977] to study the radial and latitudinal gradients and temporal variations of 
anomalous O. Quiet-time data from 1977 to 1985 are utilized. 

2. Observations. The general features of the 1 1-year cosmic-ray modulation cycle may be exam- 
ined by using the counting rate of particles which penetrate the high-energy telescope on Vl as 
shown in Fig. 1 (curve ’P". mainly protons >75 MeV). A change in this rate by a factor —2 is 
observed between 1977 and the minimum in 1981-82. This change is mainly due to solar 
modulation; however, some gradient effects are present in the data as well since VI is moving 
outward from 1 to — 22 AU during this time period. Figure 1 also shows a similar plot for 40- 
!06MeV/nuc C+O nuclei for 24 quiet-time intervals between 1977-85 (curve ’C+O"). Here the 
change in intensity between 1977 and the minimum in 1981 is a factor —5. Finally, in Fig. 1 
the temporal variations for 5.6-17.2 MeV/nuc anomalous O are shown for the same 24 quiet- 
ume intervals (curve “O’). All three components she w a similar pattern of variations but for the 
anomalous O the overall intensity variation is larger tnan the others, approaching a factor — 100. 

In order to minimize contamination by sols>r and interplanetary particles and to help 
separate temporal and radial variations, six quiet-time intervals labeled A-F in Fig. 1 were select- 
ed by setting limits to the maximum low-energy helium flux, in a manner similar to Mat 
described in Cummings et al. (1984 ]. The anomalous O spectrum is derived from the observed 
total O spectrum for a particular time interval by subtracting both a low-energy solar or inter- 
planetary component and the high-energy galactic cosmic-ray component [Cummings et al.. 
1985], 

In Fig. 2a. b. c wc show the differential intensities of anomalous O nuclei in three separate 
energy intervals for each of the six qui-t-time intervals at each of the spacecraft as they moved 
outward from the sun. The solid lines are least-squares fits to the data points, except for the 
lowest energy interval for periods E and F when Vl is at —24 degrees north heliographic lati- 
tude. (Vl left the ecliptic plane after encounter with Saturn in 1980.) In these two cases the Vl 
points, connected by the dashed lines, are significantly above the straight lines (representing a 
constant radial gradient) connecting the V2 and P10 points. We attribute the deviation to a lati- 
tude gradient at 'ow energies having a weighted average value of 3.0 ± 1.0%/degree for the com- 
bined two periods E and F. 

In Fig. 3 we show the calculated radial gradient (slope of solid lines in Fig. 2) for each of 
the three energy intervals. In the lowest energy interval (7. 1 - !0.6MeV/nuc) the spectrum of 
anomalous O is changing with time [Cummings et al.. 1985], which may contribute to the ob- 
served variation in the radial gradient. Above 10.6 MeV/nuc no such energy spectral changes are 
found and the resulting gradient is remarkably constant. In the 10.6-17.1 MeV/nuc interval the 
average value of the gradient is 10.7 ±0.6%/AU, somewhat lower than the 15 ± 3%/AU found by 
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Fig. 1. Counting rates of three particle types from Voyager 1 telescopes. P labels the penetrating 
rate from the high-energy telescope (mainly protons > 75MeV), C+O labels the rate of 40- 106 C+O, 
and O labels the rate of 5.6- 17.2 MeV/nuc anomalous oxygen. The approximate times of the solar 
magnetic field reversal in the northern (N) and southern (S) polar regions are indicated by the vertical 
bars [Webb et al.. 1 984 j. The horizontal bars represent quiet-time intervals as discussed in the text * 


0 



Radial Distance (AU) Racial Distance ( AU ) Radial Distance (AU) 


Fig 2. Intensity of anomalous O versus heliocentric radial distance in three energy intervals from 
the VI, V2, and P10 cosmic-ray instruments. The labels A-F refer to the six quiet-time intervals 
chosen for analysis. The lines are described in the text 
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Fig. 3. Measured radial gradients of anomalous O 
versus time for the three energy intervals of Fig. 2. 


Webber et al. [ 1 98 1 ] for Thr period 1972-79. 

It is clear from Figs. 2 and 3 that the radial 
gradient for all three energy intervals has not 
changed appreciably with time or radial dis- 
tance. maintaining a value between ~ 10 - 
15%/AU despite the fact that the intensity 
has varied by more than a factor of 100 . 

3. Discussion. We believe that the constancy 
of the gradient with respect to time both be- 
fore and after the solar magnetic field rever- 
sal in 1980, and during a period u-hen the in- 
tensity change was a factor — 100 is an im- 
portant clue and constraint on the solar 
modulation process. In conventional modu- 
lation theory, if the effects of p ' ,, cle stream- 
ing in the interplanetan .uedium can be 
neglected, then the radial gradient, G P and 
the radial diffusion coefficient, K p are related 
by G r -CV/K r . The average solar wind velo- 
city, V, and K r change only slightly over a 
solar cycle (see Hedgecoek 11974] and Feld- 
man el al [1979]). Also the Compton- 

Getting coefficient C, which is related to the spectral shape, does not change appreciably since the 
spectral shape remains similar above — 10 MeV/nuc. Thus, if the gradient can be described by 
these local parameters, it would be expected to be constant as observed. 

In this simple conventional modulation model in which K r is independent of radial dis- 
tance r. the particle intensity j at r is given by j-j b explCV(r-r b )/K r ], where j b is the intensity at 
the modulation boundary. Since C. V, and K r are not changing appreciably with time, and as- 
sumes, j b is constant, the large modulation at a given position would require a time variation in 
the boundary location r b . The magnitude of the required change in the boundary distance can 
be estimated from Fig. 2b. The dotted line is an extrapolation of the intensity in the 10.6- 
1 7. 1 MeV/nuc energy range during the time of solar minimum to the position of P10 (27 AU) 
during the solar maximum period of interval D. The implied intensity change from period A to 
D at 27 AU is a factor of 1 36 ± 52, indicating that most of the modulation during period D is oc- 
curring beyond 27 AU. Using the average radial gradient for this energy interval we find that the 
required change in the boundary distance r b is 46 ±4 AU. A similar boundary shift has been sug- 
gested by Evenson et al [1979] to explain electron observations. 

Such a boundary shift would also produce changes in the intensity of other particle 
species. For example, the penetrating particle rate (P) in Fig. 1 shows a variation from solar 
minimum to solar maximum of a factor —2. The median rigidity of the particles dominating 
this rate is — 1 . 8 GV, the sam* rigidity as 7 MeV/nuc anomalous O if O is singly charged (as ex- 
pected if they are freshly-ionized neutrals [Fisk et al., 1974]). The spectral shapes of the 
anomalous O (above lOMeV/nuc) and high-energy protons are similar (spectral index —-2) im- 
plying they have similar values of C. Since K r oc(JRR), where fi is the particle velocity and R is 
the rigidity, the expected gradient for the P rate can be scaled from the gradient of anomalous O 
by G/P)“G,(OKPo/Pp) ■ — 10 (. 12/.89) -1.5%/AU, a value in approximate agreement with ob- 
servation: by others (sec McKtbben et al [1982] and Lockwood and Webber [1984]). The inten- 
sities at two different times are related by j i/j 2 *cxp(G r ( 8 rb)). Therefore, a boundary change of 
46AU and a radial gradient of 1.5-3%/AU would result in j 1 /j 2 “ 2 - 4 for the penetrating parti- 
cles. similar to what we observe. 

An alternative way to accomplish the same modulation without changing the boundary 
distance is to decrease K p and thus increase the gradient, in a localized shell of turbulence in the 
outer heliosphere, as suggested by Burlaga et al [1984], The observed change in modulation 
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would require, for example, an increase in the anomalous O gradient from 11%/AU to 44%/AU 
in a shell 1 5 AU thick. 

It is also possible that variations in the source of the anomalous O may contribute to the 
observed solar cycle variation. We note that Jokipii [19851, using a model of the anomalous 
component in which the particles are accelerated at the polar termination shock and drift to the 
solar equatorial regions, has calculated spectra of anomalous O that depend on the polarity of the 
magnetic held and which are in reasonable agreement with observations of changes in the spectra 
associated with the magnetic field reversal (see Cummings el al. (19851). However, we find no 
evidence of a change in the gra'iient of the anomalous O at the time of the solar magnetic field 
reversal as would presumably be etyected from such a model. 
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Previously, only one measurement of electrons at a location remote from earth has been 
reported [Eraker, 19821. That measurement was made with a cosmic ray instrument on board 
the Pioneer 10 spacecraft during the period from 1972 - 1980, when Pioneer proceeded outward 
from 1 to 21.5 AU. Electrons of 1.75 - 25 MeV were observed. This energy range appears to be 
dominated by Jovian electrons, and no significant temporal or radial effects due tc galactic elec- 
trons were observed. 

In this paper we report on the flux of electrons with energy from — 10 - 180 MeV meas- 
ured with the electron telescope on the Voyager 1 and 2 spacecraft [Stone et ai. 1977) in the 
heliocentric radial range 1 - 22 AU between 1977 and 1984. Jovian electrons were clearly ob- 
servable between 1978 and 1983 (radial range 2 - 12 AU) at energies below — 50 MeV [Christon 
et al.. 19851 Above — 50 MeV the election intensity exhibited temporal variations generally re- 
lated to the 1 1 year modulation of protons >75 MeV. The overall magnitude of the electron 
intensity changes between the maxim”m intensity observed in 1977 and the minimum intensity 
in 1981 was a factor — 2, also comparable to that observed for > 75 MeV protons. By early 
1985 the electron intensity had apparently recovered to the level observed in 1977 whereas the 
proton intensity was still about 20% lower. A detailed interpretation of these electron variations 
in all energy channels depends on an accurate subtraction of background induced by energetic 
protons of a few 100 MeV. This subtraction is facilitated by calibration results at several ener- 
gies. Further results on these temporal variations and limits on possible interplanetary gradients 
will be reported. 
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